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Abstract 



We study the characteristics of correlation between particles in jets produced in 

heavy-ion collisions. In the framework of parton recombination we calculate the r\ and 

cf) distributions of a pion associated with a trigger particle. The origin of the pedestal in 

At] is related to the longitudinal expansion of the thermal partons that are enhanced by 

the energy loss of hard partons traversing the bulk medium. The peaks in Ar] and A<j) 

are related to the same angular spread of the shower partons in a jet cone. No artificial 

J> ' short- or long-range correlations are put in by hand. A large part of the correlation 

between hadrons in jets is due to the correlation among the shower partons arising 

from momentum conservation. Recombination between thermal and shower partons 

lO ' dominates the correlation characterisitics in the intermediate pt region. 

O ' 

in 

O ■ PACS numbers: 25.75.Gz, 25.75.-q 

1 Introduction 

The discovery of the medium effects on jets produced at the relativistic heavy-ion collider 
(RHIC) has contributed greatly to the understanding of the physics underlying hard partons 
traversing a hot and dense quark-gluon systemfl, 2]. Subsequently, a large number of exper- 
imental investigations have revealed details of the properties of jets, both on the near side 
and on the away side, in heavy- ion collisions (HIC) [3]- [8]. Particle correlations within jets 
are shown to have distinctive properties in the azimuthal angles and in pseudorapidities 
r] of the associated particles. Of particular noteworthiness is the peak in Ar] between the rj 
values of the trigger and the associated particles on the near side: it sits above a flat plateau 
or a pedestal [3]. No pedestal is found in A0, although a peak in that variable is prominent. 
It is our aim in this paper to study the peaks and pedestal in A0 and Ar? and to discuss the 
physical origins of these jet characteristics. 

Our approach to the problem of correlation between particles in jets will be based on 
parton recombination that has been shown to be successful in reproducing the single-particle 
distributions in Au+Au collisions [9], and in d+Au collisions [10]. It has also been applied to 
the study of dihadron correlations [11, 12]. However, since those investigations were carried 
out in the framework of the one-dimensional (ID) formulation of the recombination model, 
a generalization to account for the 3D geometry of the problem is necessary if we are to 
describe the features in the and r] variables. Once we consider the 3D aspects of jets 
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the angular distribution of the shower partons in the jet cone becomes an important new 
feature. Our treatment will be phenomenological, since no results from perturbative QCD 
considerations can reliably be applied to the intermediate pt region, which is where the 
data on correlations are taken. Furthermore, the pedestal phenomenon suggests that the jet 
interacts with the environment in a way that is sensitive to the hydrodynamical expansion 
of the medium, since there is longitudinal expansion, but not azimuthal expansion. It turns 
out that we can relate the pedestal to the local enhancement of the thermal bath due to the 
energy loss of the hard parton to the medium. 

Our approach to the problem, though not based on first principles, provides the first 
theoretical interpretation of the pedestal and peak structure of the jets produced at RHIC 
in quantitative terms. The roles that shower and thermal partons play in the hadronization 
process are found to be crucial and cannot be reinterpreted in any other model that we are 
aware of. 

In Sec. 2 we describe how the 3D properties of jets are incorporated in the recombination 
formalism. In Sec. 3 the pedestal phenomenon is investigated, followed by the study of the 
peaks in Arj and A0 distributions in Sec. 4. The conclusion is given in the last section. 

2 Two-Particle Distribution 

In the framework of parton recombination formulated in ID for particle production in HIC 
[9]- [12], let us start by writing down the single-and two-particle distributions 

— = ~2 [— — ^2(91, 92)^1(91,92, Pi) (1) 
pxdpx pf J qi q 2 



! ' — } F 4 (q l ,q 2 ,q 3 ,q 4 )R 7T1 (q 1 ,q 3 ,p 1 )R 7T2 (q 2 ,q A ,p 2 ) , (2) 



Pxp 2 dp l dp 2 p\p\ J \^ qj 



where the recombination function (RF) for a pion is [13] 

^(9i,92,p) = ^f- + -- 1 ) • ( 3 ) 
p 2 \p p J 

The two- and four-parton distributions, F 2 and F 4 , can be written in terms of their compo- 
nents as 

F 2 {1,2) = (TT + TS + SS) 12 , (4) 

F 4 (l, 2, 3, 4) = (TT + TS + SS) 13 (TT + TS + SS) 24 . (5) 
The thermal parton distribution is 

T(q) = Cqe^ T , (6) 
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and the shower partons in a jet have the form 



S(q)=tY,f dkk M k ) s 'iW k )> (7) 

i J 

where fi(k) is the distribution of hard parton % in a HIC, and £ is the average fraction of 
hard partons that emerge from the thermal medium to hadronize. S{ is the shower-parton 
distribution (SPD) of parton j in a shower initiated by a hard parton i; their parameteriza- 
tions for various % and j are given in [14]. For the sake of clarity we have omitted the labels 
for parton species for T and S in the above equations and in what follows, except when the 
identity of a parton is needed for emphasis. 

It is important to recognize that Eq. (7) should not be applied to every S that appears 
in Eqs. (4) and (5) because all the shower partons in F 2 and F 4 are created in the same jet, 
so all the SS ■ ■ ■ terms share the same hard parton with momentum k. Thus for two shower 
partons (SS), their joint distribution is 



(ss)( qi , q2 ) = ce/ dkk m hi (j^h)} , (8) 



where the quantity in the curly brackets has the form in terms of momentum fractions X\ 
and X2 



{st M , s{ (^_)} = i [s!( Xl) s( (^-) + sj (^-) 8f M 



(9) 



This is done to guarantee momentum conservation x±+x 2 < 1 and to symmetrize the order of 
emission of the two partons. Not only should Eq. (8) apply to two adjacent shower partons, 
like in (<S<S)i3 and (<S<S) 2 4, which recombine to form the pions at p\ and p 2 , respectively, 
but also to two shower partons like in (TiS>) 13 (TiS>) 24 , which do not recombine with each 
other, but separately with thermal partons. It is the structure in Eq. (8) that gives rise to 
correlation between two shower partons in a jet, and thereby endows the detected hadrons 
at high pt with correlation. However, that is sufficient primarily in the case of considering 
only the correlation among momentum magnitudes in ID. 

To generalize our consideration to 3D, let us focus our attention first on the (TS) com- 
ponents in Eq. (5) that give the most important contribution to the trigger and associated 
particles in central collision. To be specific let us further designate the trigger particle to be 
a 7r + and the associated particle to be a n~ . The 3D expression for Eq. (2) in that case is 



dN ^- = 1 ,y f d]±m f (f\d 3 r d3qj 

Pxdp x dr) tr d§ tr p 2 dp 2 drid<\> (P1P2) 2 i' J L = i 3 Qj 



xSf (fi, cfi) T d (f 3 , q 3 ) R { ^1 (fi, f 3 , r pi ; gi, (f 3 ,Pi) 

xSf (f 2 , q 2 ) T n (f 4 , q 4 ) R { *1 (f 2 , f 4 , f p2 ; q 2 , q 4 ,p 2 ) (10) 

where is the 3D version of the RF and it includes the spatial coordinates. Since re- 
combining partons must not only have collinear momenta but also have overlapping wave 
functions, there are many narrow Gaussian distributions in R^ that reduce the phase space 
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of integration. Let us first restrict our attention to the x-z plane, i.e. the plane containing 
the hard parton momentum k and the longitudinal direction z, and assume that all q} vectors 
are in that plane so that we may write 

ft = (qi, 6i), Pl = (Pl, Otrig), P2 = (P2, 9) (H) 

and 

7fc = In cot Vtrig = hi cot 7] = In cot ^. (12) 

We further assume that qi is in the direction of k and that the trigger momentum pi is also 
along fc. This is a simplification that does not compromise the angular correlation between 
the trigger and associated particles, which we shall study in detail. Since requires the 
wave functions of the recombining partons to overlap, the shower and thermal partons at T\ 
and r 3 should be nearby. These assumptions allow us to reduce the trigger momentum part 
of Eq. (10) to the usual ID formulation 
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j3 j3 d 3 q 1 d 3 q 3 a-dt- -\r>{S)r- - - - - - \ 

d rid r 3 S i (r^qJT (r 3 , q 3 ) R y J (r u r 3 ,r pi ; q u q^p^ 

qi q3 



1 



qi q 3 \k 



R$(qi,q S ,Pi), (13) 

-8l=9 tr ig 



where R^ is the usual RF in ID, given in Eq. (3). The superscript d on T is omitted on 
the RHS since the flavor dependence of the thermal partons is negligible and ignored. We 
shall assume that at midrapidity, —0.7 < r\ trig < +0.7, the thermal distribution does not 
contain essential dependence on i] 3 so that the condition 63 = 61= 6 tr i g does not lead to any 
restriction that would cause T(g 3 ,# 3 ) to deviate from the usual parameterization given in 
Eq. (6). 

In the case of the associated particle we must recognize that the shower parton q 2 may 

— * — * 

be emitted at an angle ip relative to the hard parton momentum k. Since k forms an angle 
61 with the z axis (after identifying the direction q\ with k), we have 

i> = 62-61- (14) 

The recombination of the shower parton at q\ with a thermal parton at q 4 requires that 
they are not only collinear, but also overlapping spatially. There are therefore 5-functions, 
5(6 2 — 6 A )5 ( 02 + 04 — that force the momentum p 2 of the associated particle to be in the 
same direction as q 2 and q 4 . The corresponding 3D integrations result in 
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j3 j3 d 3 q2 d 3 q 4 , _ u ( 3 ) _ _ _ _ _ _ 

d 6 r 2 d 6 r A S t (r 2 ,g 2 )T (r 4 , q 4 ) R^l (r 2 , r 4 , r P2 ; q 2 , q 4 ,P2) 

<?2 94 



/ 



<?2 <?4 \k-qi J 



R { - ] -{q2,q^P2) ■ (15) 

6> 4 =02=01+^ 



This equation is very similar to Eq. (13) but with one crucial difference; that is, the shower 
parton distribution now depends on ip. In Eq. (13) no angular dependence of the shower 



4 



parton gi is assumed for simplicity, since it is only the difference between q\ and that 
matters. That angular difference, vp, made explicit in Eq. (14), is translated into the angular 
difference between pi and f>2- We have tacitly assumed that the thermal partons have local 
angular spread that allows a q% to recombine with q\ to form pi, and similarly q± with q 2 to 
form p 2 , in such a way that pi and p*2 are in the directions of q\ and 52, respectively. 

The dependence of a SPD on ip cannot be determined from pQCD, since the values of 
p T concerned is in the 1 < p T < 4 GeV/c range. We shall assume that the SPD in Eq. (15) 
can be written in the following factorizable form 



12 



-A 
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\k-qi J \k-q! 
where the dependence on ip has a Gaussian distribution 

G ^ x) = eXP [-2^{x)_ ' 
whose half-width depends on the momentum fraction x of the shower parton 



(Toil 



X) 



(16) 



(17) 



(18) 



with <jq being a parameter that is adjustable. Sf on the RHS of Eq. (16) is to be sym- 
metrized with Sf in Eq. (13) without further complication from angular consideration. The 
momentum fraction relevant for a(x) is independent of q\. Once we obtain Eq. (16) 

where G(ip,q2/k) stands as a modifying factor, it is reasonable to liberate that factor from 
the way in which it is derived, and allow it to assume a 3D property so that when we later 
consider the A0 behavior in the transverse plane the same factor G(ip,q2/k) applies to the 
azimuthal angle between pi and p 2 with ip replaced by A<f>. Using Eqs. (12) and (14) we 
obtain with 9 = 9 2 



tan - = tan^— = g{rj : 771) = - 



v _ 



-VI 



-r)-m 



(19) 



where rji is to be identified with rjtrig- 

We now can write Eq. (10) with the help of Eqs. (3), (13) and (15) in the simpler form 
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p 1 dp 1 d,7] 1 d(j) 1 p2dp2drid(f) (P1P2) 3 



^ / dkkfi(k) J dq 1 J dq 2 9(k - Ql - q 2 ) 

T(p 2 - q 2 ,v) 



3,3 



Qi 



Si I — I , Sj 




tp=2tan 1 g(f],rn) 



(20) 



where i is summed over all hard parton species. If i is a valence quark of tt + or n~ , the 
corresponding SPD is K, while if % is a sea quark, the corresponding SPD is L in the notation 
of [14], where K = S^ al+sea and L = S* ea . The integral over k will, in practice, be from 3 to 
30 GeV/c, and will be integrated from to p«, although the formalism is not reliable for 
any transverse momentum less than 1 GeV/c. There are other terms of F4 contained in Eq. 
(5) that can be written out as in Eq. (20), but will not be detailed here. 
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3 Pedestal and Background Subtraction 



The experimental procedure of making background subtraction involves multiple consider- 
ations, and the result is that there remains a residual pedestal in A77 in central collisions, 
but no pedestal in A<f>, where A77 and A0 are differences in 77 and between the associated 
particle and the trigger. In [11] the dihadron correlation is calculated with the assumption 
that the factorizable part of the two-particle distribution corresponds to the background. 
That leaves only the (TS + SS)i%(TS + SS)2i terms of F 4 to consider. We now find that 
in order to understand the origin of the pedestal what constitutes the background must be 
reconsidered. 

The first point that we want to address is the effect of energy loss of hard partons travers- 
ing the hot medium. The parameter £ in Eqs. (7) and (8), being only 0.07 as determined in 
[9] for central Au+Au collisions, implies that a large fraction of the hard partons produced 
are absorbed or attenuated by the bulk medium. Those that emerge to hadronize outside 
must have been created near the surface. In the short distance that such a parton travels in 
the medium it must on average lose some energy and locally enhance the thermal motion of 
the partons in the environment. Those enhanced thermal partons should have an effective 
inverse slope T that is slightly higher (by AT) than that of the usual thermal partons not 
influenced by the passage of any hard partons. The latter corresponds to what is measured 
for the low-p T pions in single-particle distribution, which for p T < 2 GeV/c has led to the 
determination of T = 0.317 GeV/c in [9]. If AT is small, as we expect it to be, it does 
not make too much difference in the calculation of the single-particle distribution in the 
intermediate pt region. However, in principle, in thermal-shower recombination the shower 
partons recombine with the enhanced thermal partons that are in the immediate vicinity 
of the hard parton that creates the shower. To distinguish the two values of the inverse 
slope, let us use T to denote the one for the enhanced thermal parton, since that is the only 
physically relevant one for TS recombination, i.e., it is what appears in Eq. (6) for T. Let 
us use T to denote the inverse slope for the thermal medium undisturbed by jet quenching, 
i.e., T = 0.317 GeV/c. We define 

AT = T - T , (21) 

which is a parameter to be determined phenomenologically in this paper, although it can 
presumably be determined directly by dedicated experimentation. We shall relate AT to the 
pedestal. The physics involved is sensible, since the pedestal is seen only for central collisions 
where TS recombination is important. In peripheral collisions thermal partons play a minor 
role and there are few of them to be enhanced. 

We now give an argument why the pedestal is seen in Arj, but not in A0. When a hard 
parton is first scattered or created at a large angle relative to the incident beam direction, 
the bulk medium is in a highly compressed state. As it expands longitudinally, the region of 
enhanced thermal partons expands with the whole system, although limited to the neighbor- 
hood of where the hard parton traverses the thin layer close to the surface. Thus a section 
at midrapidity can have the enhanced T. In the transverse plane, on the other hand, the 
expansion is in the radial direction only; that is, there is no mixing across different azimuthal 
sections. The subtraction scheme carried out by STAR [3] defines the yield dN/dAip to be 
zero at |A0| = 1 after subtraction. Since there is no expansion in the <fi direction, we can 
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take this subtraction scheme to imply that all thermal partons that are enhanced stay inside 
the \Acf>\ < 1 region, and that for \A(f>\ > 1 there is no thermal-thermal recombination that is 
not in the background. There may still be an enhancement of the yield inside the \A(f>\ < 1 
region that plays a role similar to the pedestal in Arj, but it cannot extend beyond \A<j>\ = 1 
by virtue of the subtraction scheme. 

The above description cannot be demonstrated by a transport model, since we do not 
have a Monte Carlo code for the evolution process. However, the physical content of our 
discussion will be embodied in our quantitative formulation below. The recombination model 
actually provides a more transparent description of the hadronization process than what a 
code without reliable equations can offer. 

Let us first consider the particles that are associated with a trigger. Since in Au+Au 
collisions the TS component is dominant in the intermediate pt region [9], which includes 
the trigger window 4 < < 6 GeV/c, we select the (TiS) 13 term in Eq. (5) for our discussion 
below, although both terms in (TS + SS) 13 are included in our calculation. Thus the parts 
of F 4 that contribute to the trigger and its associated particles are 

Fl r+as = (TS) l3 (TT + TS + SS) 24 . (22) 

Among the three terms on the right side it is clear that (TS + SS) 24 are directly related 
to the jet and contribute to the peak in the associated particle distribution (APD), while 
(TT) 24 does not involve the shower partons but can nevertheless contribute to the APD 
outside the peak. It is therefore our prime candidate for the pedestal. The absence of an 
obvious pedestal in the APD in A<p is a consequence of the background subtraction scheme 
in the analysis of the experimental data, where the A0 distribution at \A(f>\ = 1 is identified 
to be the background. In our description of the APD the background corresponds to 

Ft 9 = (TS) 13 (T T ) 24 , (23) 

where %(q) is the thermal parton distribution in the absence of any hard partons, i.e., 

%(q) = Cqe~ qlTa . (24) 

The meaning of To has already been discussed in connection with Eq. (21). We have tacitly 
assumed that the normalization factor C is unchanged from that in Eq. (6). The experimental 
subtraction scheme implies that 

(TS) n (TT) 24 = (TS) 13 (T T ) 2A (25) 

at \A(f>\ = 1. We now define the APD with background subtracted to be generated by 

Ff p = (TS) 13 (TT-T T + TS + SS) 24 , (26) 

in which (TT — T T ) 24 vanishes at \A(f>\ = 1, but need not be zero at \A(f>\ < 1. Equation 
(26) is our basic input in the calculation of the APD. 

Let us now write our two-particle distribution in a form suitable for application of the 
experimental cuts. In Eq. (20) we have a differential distribution in dp\dr\xd$\dp 2 dr}d$ where 
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rji and 0i denote the variables of the trigger, and rj and those of the associated particle. 
With the definition 

A77 = 7/ - 771, A0 = 0-0x , (27) 

we change the differentiated variables to dpidp2di] 1 dAi]dA(f), with <pi being used as a free 
reference point in the assemblage of the A0 distribution; the trigger rji is to be integrated 
over the trigger window —0.7 < 771 < +0.7. We then have 

dN 

iV(l,2) ^ 



dpidp2dr]idAr]dA(f) 
4 dq 



= /(n — | F 4 P ( ( lu ( l2,q3,q4)R7r 1 ( y ( li, ( l3,Pi)Rn 2 ( y q2,qA,P2) ■ (28) 

P1P2 •> y i=1 ?j y 

The experimental APD involves integrations over the trigger window in r]\ and p\ and the 
window of the associated particle in p 2 . It is 

dN AP g:dp 2 f?d Pl f jl 7 d Vl N(l,2) 

__ _ __ ; ^yj 



where iV(l) is the one-particle trigger distribution that receives contribution from the (T<S) 13 
component only, i.e. 

d N tr 1 - c ^ / n \ 

^ (1) = rf^ = ^^/^ /l( ^ gi ^ T(pi_gi ^ l)5 Hl) • (30) 

The two-particle distribution iV(l,2) in Eq. (29) can be divided into two explicit pieces. 
For the remainder of this section we consider only the piece related to the pedestal 

JATped 

N(l 2) ped = " ' 



dp 1 dp2dr]idAr]dA(p 
= ( pi ^ 2 ^E / dkkf l {k)dq 1 dq 2 

£ Tip, - qi , Vl )Si (|) [T(q 2 , A V )T(p 2 - q 2 , A^)iJ(A0) 

-%{q 2 )T Q {p 2 -q 2 )} , (31) 

where we have extracted the <f> dependence of the associated particle in the form of a factor 
H(A(f)), which satisfies some constraints to be specified below. It should be noted that Eq. 
(31) is factorizable, since the two parts of (T5) 13 (TT — T T ) 24 are independent of each 
other. The (TS)is part contributes to the trigger distribution given in Eq. (30). It will be 
cancelled in the ratio defined in Eq. (29). 

Using Eq. (24), the integration of the last term in Eq. (31) over q 2 can readily be carried 
out, giving 

1 fP* , , C 



rp2 (J* 

J o dq 2 T {q 2 )T {p 2 -q 2 ) = —exp{-p 2 /T ) , (32) 



which, upon further integration over the window (p a ,Pb) of the associated particle, results in 

dN , . , ■> - i tut rn \ (33) 



rpb 

J Pa Pldp2 



l (rr \ 2 h( Va Pb 



where 



h(x, y) = (1 + x)e x - (1 + y)e 



(34) 



Similar result follows for the TT term in Eq. (31). 

Next, we consider the Arj and A(f) dependences of the pedestal in Eq. (31). As we have 
discussed qualitatively in the beginning of this section, the longitudinal expansion allows the 
enhanced thermal partons to extend over a wide range in 77, but radial expansion does not 
increase the range in 0. Putting these properties into quantitative terms, we give T(q 2 ) no 
essential dependence on A77 in Eq. (31), but a Gaussian dependence on A0 so that 



H(A<p) = c exp ( -A(p 2 /2a, 



(35) 

The normalization c is to be determined by the condition that when A0 is integrated over 
the experimental window (—0.5, +0.5) we obtain 



/°' 5 dA<j)H(A<p) = 1 
J-0.5 



(36) 



In this way we can relate Eq. (31) to the pedestal observed in the experiment without explicit 
factors dependent on the experimental windows. The half-width is adjusted so that the 
quantity inside the square brackets in Eq. (31) vanishes at |A0| = 1, as required by the 
subtraction scheme, Eq. (25). What we do here is to fix all the extra free parameters of 
the problem by the experimental features of the data. Since our aim is to reproduce the 
observed characteristics of the data, which are presented with specific cuts, it is impossible 
to do so without incorporating those cuts. However, it does not imply that we are merely 
fitting the data with free parameters. We shall perform several multi-dimensional integrals 
to obtain in the next section the APD with the peak whose magnitude is a prediction of our 
model calculation. 

Upon integration of Eq. (29) over A0 in the acceptance window, we get for the pedestal 
part 



dN ped 

dAr] 



0,5 dNPed 

dA(j) - — 
-0.5 dArjdAcp 



C 2 
~~6 



T 2 h 



T ' T 



T 2 h 



Po_ Pb_ 

T ' T 



(37) 



where we have made use of the factorizability of Eq. (31) and the result of integration over 
P2 given in Eq. (33). This is the constant pedestal in Arj, which we relate to AT. Using 
C = 23.2 (GeV/c)- 1 , T = 0.317 GeV/c from [9], p a = 2, p b = 4 GeV/c as in the experiment 
[3], and Eq. (21) for the definition of AT, we fit the pedestal height of ~ 0.05 [3] and obtain 



AT = 15 MeV 



(38) 



Thus the inverse slope of the enhanced thermal partons is only 5% higher than that of the 
un-enhanced partons, small enough to have a negligible effect on the calculation of the single- 
particle distribution using T . However, the effect on the associated particle distribution is 
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evidently not negligible. Although Eq. (38) is a fitted result, it should be emphasized that it 
is in the specific model of parton recombination that the pedestal phenomenon is interpreted 
as a consequence of the enhancement of the thermal energy of the soft partons due to energy 
loss of the hard parton. 

Equation (37) is the projection of dN AP / dAr/dAcj) onto Ar/. We now project it onto A<p 
by integrating it over A 77 from —1 to +1 in accordance to the experimental cut [3] and obtain 
from Eq. (31) for the pedestal part 



dN ped 
dA<f) 



J dAr] 



dN ped 

dAr]dA(f) 



3 



T 2 h 



Pa Pb 

T ' T 



H{A(f)) - 1%h 



Pa Pb_ 

T ' T 



(39) 



This is not a flat pedestal as in the case of A77, but is nevertheless a remnant of the effect of 
enhanced thermal partons sitting under a peak from the (TS + SS) 2 a part of F± p that we 
have not yet calculated. Although it is not visible as a plateau, its origin is the same as that 
which gives rise to the pedestal in A77. It is the background subtraction in A<fr that forces 
Eq. (39) to vanish at |A0| = 1, which is achieved by our choice of 0$ in Eq. (35) its value is 



0> 



1.2 



(40) 



With this choice of 0$ Eq. (39) is entirely specified numerically. The A0 dependence of 
dN ped /dA(j) is a broad mount between A0 = ±1, which we shall exhibit later together with 
the complete dN/dA(j). This mount is the projection of a ridge in the 3D display of APD in 
(A77, A<p) onto the A0 subspace. 

This completes our discussion of the (TT — T0T0) contribution to F± p in Eq. (26) and 
then to Eqs. (28) and (29). We now proceed to the last two terms of F± p '. 



4 Peaks in Ar] and A0 

We now return to Eq. (26) and consider the (TS + SS) 24 terms in F 4 AP that we have 
put aside. The TS component has already been described in detail in Eq. (20). The 
SS component is less important except in peripheral collisions. It can be included by the 
replacement of the {Sj, Sf j T term in Eq. (20) by a {s{, S{ , Sf j term with appropriate 
symmetrization. However, since the recombination of two shower partons reproduces the 
fragmentation function in accordance to 

xD?(x) = I ^^{si( Xl ),S{ (-^-)}R n ( Xl ,x 2 ,x) , (41) 

J X\ X2 I \ 1 — X\ / J 

from which the SPDs are derived in the first place, we may replace a SSR term by a xD 
term. For notational simplicity we leave the (iSiS>) 2 4 terms out in our description below, but 
numerically include their recombination in the calculated result. 

Let us recall that starting from the general formula Eq. (10) we first focussed on the 
variables in the plane containing k and the longitudinal direction z and related the relevant 
angles to pseudorapidities in Eq. (12). We then defined the angle ip between the shower 
parton q 2 and the hard parton k to be as given in Eq. (14). The angular dependence of the 
SPD is then expressed by a Gaussian distribution G(ip,x) in Eq. (17), put in a factorized 
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form in Eq. (16). We may identify G(if),x) as the x-z projection of a general Gaussian 
distribution that describes the dependence on the angle between pi and p 2 in 3D geometry. 
In the small width approximation, it can be shown that the corresponding projection onto 
the transverse, x-y, plane, gives a Gaussian form, G(A<f>,x), for the A(f) dependence. 

— * 

Returning to Eq. (20) which is written for the plane containing k and z, we now write 
the double differential in both Arj and A0, as for the pedestal term, but now for the peak 
term arising from the (T<S)i 3 (T<S) 2 4 component 



dN peak_ 



—^CYl J dkk fi( k ) J d qi J dq 2 9(k -qy- q 2 ) 



dpidp 2 dr]idAridA(j) (P1P2)' 

£T( Pl -< h) {s!m.s![-^-)\T {l ,- <l2) 



n 

k 




i/>=tan 1 g(r],r]i) V k 



bG ( A0, , (42) 



where b is a numerical normalization factor such that when Eq. (42) is integrated over A0 
from —0.5 to +0.5 (the experimental window for projection to the Arj dependence) one gets 



b f ' dA(f)G(A(f),x = 0) = 1 (43) 

J-0.5 

in the small x approximation. 

Identifying Eq. (42) as N(l,2) peak , we can substitute it into the numerator of Eq. (29) 
and obtain dN peak / dAr]dA(p. There is only one free parameters to adjust; it is cr in Eq. 
(17). The projections of the double differential distribution to Arj and A0 separately can be 
compared with the data when the pedestal component is added. That is 



dN AP r°> 5 



dAr] 



/U,0 
dA(f) 
-0.5 



dN ped dN peak 



dAr]dA(p dAr]dA(j) 



(44) 



dA<p 



J dAr] 



dN ped dN peak 
+ 



dN AP 

(45) 



dAr]dA(p dAr]dA(j) 



The pedestal contributions to the above integrals are given by Eqs. (37) and (39), respec- 
tively. It should be noted that the latter has a A0 dependence described by H(A<f>), which is 
a Gaussian with a half-width given by Eq. (40). We shall find below that the half-width of 
the peak term, controlled by o"o, is much smaller than a^, so the last terms in Eqs. (44) and 
(45) dominate the peak structure of the APD. The term dN ped / dAr]dA<p gives rise to a ridge, 
whose projections are the flat pedestal in dN AP /dAr], and a broad mount in dN AP /dA(p. 

After we put all the pieces together, using Eq. (29) as the important link between our 
calculation and the measurable quantities, the results from Eqs. (44) and (45) can be com- 
pared to data with a adjusted to fit the peak width. Figure 1 shows the APD in Arj. The 
solid line is the result of our calculation when we set 

a = 0.22 . (46) 
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Evidently, it reproduces the data [3] satisfactorily, both in the peak structure and the 
pedestal. It should be emphasized that although the width and the pedestal height are 
adjusted to fit by our choice of o"o and AT, the height of the peak is a consequence of our 
very complicated calculation involving multiple integrals over many terms based on the re- 
combination model. It is by no means trivial that the data can be fitted so well. In point of 
fact, we have only calculated 7r + 7r~ production, whereas the data are for all charged parti- 
cles. Thus the precise values of the parameters are not as significant as the overall situation 
where the pedestal and peak structure of the A77 distribution of the data can be reproduced 
by our description of dihadron correlation. 

Turning now to the APD in A<f>, we have no free parameters to adjust, since the double 
differentials in the square brackets in Eqs. (44) and (45) are identical. Upon integration over 
At] we obtain the solid line in Fig. 2, which compares well with the data [3]. The dashed 
line indicates the pedestal contribution from Eq. (39). It plays the role of the flat pedestal 
in Fig. 1, but here it vanishes at | A(f>\ = 1 because of the subtraction scheme. Since the peak 
structure in A0 is much narrower than the broad mount of the pedestal (due to a <C a^), 
we have applied the vanishing condition at \A(f)\ = 1 only to the dN ped / dA<f> component, 
knowing that the dN peak /dA<p component is negligible at the wings of the peak. 

Figures 1 and 2 represent our main results on dihadron correlation in central Au+Au 
collisions at -/s = 200 GeV, when the trigger particle is kept within the range 4 < p trig < 6 
GeV/c and the associated particle in the range 2 < p asS oc < 4 GeV/c. We can readily 
calculate the APD at higher momentum ranges, but not at lower momenta, since our model 
is not reliable for parton momentum less than 1 GeV/c. For that reason we do not calculate 
the APD for 0.15 < p asS oc < 4 GeV/c, even though data for that range are presented in [3]. 

5 Conclusion 

To summarize, we have successfully reproduced the data that show the existence of peaks 
in Arj and Acf) distributions, the former sitting above a flat pedestal, while the latter sitting 
above a broad mount. Some parameters are used to fit the data, but the essence of our work 
is not data fitting. We have demonstrated that the physics underlying the detail structure 
of the jet characteristics observed in RHIC experiments can be understood in the framework 
of parton recombination. Jets produced in heavy-ion collisions create shower partons that 
are in the environment of thermal partons, which can themselves be enhanced by the pas- 
sage of hard partons through the medium. That view is probably shared by all theoretical 
approaches to the problem. The issue then is how those partons hadronize. Different models 
treat the hadronization process differently. Our approach emphasizes the recombination of 
the thermal partons and the shower partons in the intermediate pt region. That is our way 
of accounting for the medium effects on jets. So far we have not encountered any obstacle 
in understanding the data in that way. Data fitting is only a concrete demonstration that 
the details of the jet structure can be quantitatively reproduced in our treatment. The 
hadronization formalism is thereby enriched by the determination of some features in the 
model by phenomenology. 

Our first discovery in this work is that the pedestal in the A77 distribution can be related 
to AT in the local thermal distribution. No "long-range correlation" has been put in by 
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hand. The pedestal is a consequence of the recombination of the thermal partons among 
themselves, which are only indirectly affected by hard scattering through the enhancement 
of T due to energy loss. Thus the pedestal is not a part of the jet, but cannot be present 
without a jet. This chain of successive connections involving energy loss, enhanced thermal 
partons, their hadronization by recombination, and the elevated APD in A77 is certainly very 
different from a model where an explicit correlation is put in by hand. 

Our second achievement in this work is the success in describing the peak structures in 
both At] and A0 distributions in terms of one angular distribution of the shower partons 
in a jet cone. No "short-range correlation" has been put in by hand. The only primitive 
correlation in the problem is that which exists among the shower partons in a jet. The 
properties of that correlation in the momentum magnitudes of the shower partons have 
already been described in [12]. Here we show how the correlations manifest themselves in 
the angular variables of the produced pions. 

There is one caveat in our analysis that should be noted. We have calculated only 7r + 7r~ 
production, whereas the data are on all charged hadrons. Thus numerically the parameters 
determined here are not definitive. Since particle identification is steadily being improved, 
the more appropriate arena for detailed matching of theory and experiment is when the data 
for the production of specific species become available. 

Despite the tentativeness of the value of AT given in Eq. (38) on account of the statement 
made just above, our assertion that the pedestal is related to the enhancement of T of the 
thermal partons remains unaffected. As a test of that assertion we propose that the proton 
associated with a 7r + trigger be measured in addition to the measurement of a 7r + or 7r~ 
associated particle. The p/n ratio of the associated particles in a 7r + triggered jet should be 
a good probe of the physical mechanism underlying both the pedestal and the peak. The 
recombination approach to the problem is well positioned to calculate that ratio, as it has 
for single-particle distribution [9]. We surmise that the pedestal part of the p/ir ratio would 
be higher than in the single-particle case because AT > 0. 

What we have done in this paper is influenced greatly by our intention to understand the 
pedestal and peak structure in the data of [3] . As a consequence, our analysis involves many 
integrals that correspond to the experimental cuts in the data. Having determined the origin 
of the pedestal and peaks in the Arj and A<f> distributions, we are now equipped to launch a 
study of the correlation problem independent of any triggers and related experimental cuts, 
on which more and more data are becoming available. 
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Figure Captions 

Fig. 1. Associated particle distribution in A 77 for 2 < < 4 GeV/c with trigger particle 
in 4 < p% %g < 6 GeV/c. The data from Ref. [3] are for all charged hadrons in the 
respective p T ranges. The solid line is the result of our calculation of tt~ associated 
with 7r + trigger. 

Fig. 2. Same as in Fig. 1 except that the distribution is in A0. The dashed line represents 
the pedestal effect in A0 forced to vanish at |A0| = 1 by the subtraction scheme 
discussed in Sec. 3. 
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